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m '■ Abstract 
O 

We introduce a new and well defined power counting for the effective field 
Q\ ■ theory describing nucleon-nucleon interactions. Because of the large NN scat- 

tering lengths it differs from other applications of chiral perturbation theory 
and is facilitated by introducing an unusual subtraction scheme and renormal- 
CJ ' ization group analysis. Calculation to subleading order in the expansion can 

be done analytically, and we present the results for both the ^Sq and 3 S\ — 3 Z?i 
channels. 
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Effective field theories are standard for dealing with strong interaction physics in the 
nonperturbative regime. The idea is to construct the most general Lagrangian consistent 
with the symmetries of quantum chromodynamics out of fields that create and destroy 
the relevant degrees of freedom. To have predictive power an expansion scheme must be 
found that limits the number of terms that occur in the effective Lagragian. In the case of 
the effective Lagrangian for pion self interactions the expansion is in powers of the typical 
momentum p ~ m n divided by a typical hadronic scale, of order the rho mass. At leading 
order, tree diagrams dominate with either two derivatives or one insertion of the light quark 
mass matrix at the vertex. At next order one includes both four derivative operators at 
tree level, as well as one loop diagrams with the leading operators at the vertices, and so 
forth. This perturbative expansion is known as chiral perturbation theory, and has also 
been applied to nucleon-pion interactions. However, attempts to apply chiral perturbation 
theory to nuclear physics, where nucleon-nucleon interactions are of primary interest |T|-[T6|, 
have encountered difficulties, stemming from the large scattering lengths in the ^So and 3 Si 
nucleon-nucleon scattering amplitudes The large scattering lengths generally imply 

finely tuned cancellations between graphs in the effective theory, with the result that a 
consistent expansion has so far proven elusive. 

In this letter we introduce a new method for applying effective field theory to nuclear 
physics. Our method uses dimensional regularization with a novel subtraction scheme which 
permits a consistent expansion procedure yielding analytic expressions for nucleon-nucleon 
scattering amplitudes 1 . We apply it to scattering amplitudes in both the ^So and 3 Si — 3 Di 
channels and compare with the results of the Nijmegen partial wave analysis |L7 . 



We begin by considering nucleon-nucleon scattering in the ^So channel in a nonrelativistic 
effective field theory with only nucleon fields, as is appropriate for momentum much less than 
the pion mass. The tree level amplitude is 

oo 

iA tree = -*(/i/2) 4 ~ D £ C 2 Mp 2n 

n=0 

= -z(/V2) 4 - D CV^) , (1) 

where p is the magnitude of the nucleon momentum in the center of mass frame; the total 
center of mass energy is E = p 2 /M + . . ., where M is the nucleon mass. In eq. (|l|) C2 n (/-0 
is a linear combination of coefficients of four nucleon operators with 2n spatial derivatives, 
fi is the subtraction point and D is the dimension of spacetime. In order to calculate the 
full amplitude we must sum a bubble chain of Feynman diagrams with the tree amplitude 
giving the vertices. The integral required is 



{2n) D ^ (E/2 - q° - q 2 /2M + is) {E/2 + q° - q 2 /2M + is) 



Dimensional regularization is the most convenient regularization scheme, as it preserves chiral 
symmetry and gauge symmetry; it also preserves Galilean invariance, which makes the Feynman 
integrals relatively simple to evaluate. However, physical results do not depend upon the choice of 
regulator. 



2 



= -M(MEn-ME-ie)^r(l^)^t^. (2) 

This amplitude has no pole at D — 4, but it does have a pole at D = 3 resulting from 
what would be a linear ultraviolet divergence in the integration if a momentum cutoff were 
used. The residue of the pole at D = 3 equals M(ME) n n/4ir. In the minimal subtraction 
scheme, MS, counter terms are added which subtract the poles at D — 4. Here we use a 
new subtraction procedure which we call power divergence subtraction, PDS , where the 
poles at D — 4 and D = 3 are both subtracted by counter terms. In the case of I n this 
subtraction amounts to adding the counter term contribution 

M(M£)> 

6In = ~ An(D - 3) ' (3) 

and the subtracted integral is 

I™ 8 = I n + SI n = -{ME) n (/i + ip). (4) 

The choice fi = corresponds to minimal subtraction. Using the PDS subtraction scheme 
the full amplitude, including loop effects is Q 

l + MC(p\ n)(ji + ip)/47r' 1 J 

The scattering amplitude A is related to the S'-matrix by 

S-l=e^-l=i(^)A, (6) 
where 5 is the phase shift. It is convenient to consider the quantity 

47T 47T 

P cot5 = ip+ m = - m ^-p, (7) 
which has an expansion in p 2 of the form 

pcot5 = -- + ^r p 2 + . . . , (8) 

where a is called the scattering length and tq is the effective range. If there is a state that 
is very near zero binding energy then \a\ will be very large, which is the case in the 



2 Some recent papers have criticized the applicability of dimensional regularization to NN scat- 
tering |7,il,14|]. Since C is an arbitrary polynomial in p 2 , over a range in p, eq. (6) gives the most 
general form for S consistent with unitarity. Clearly dimensional regularization can be used for NN 
scattering. Relativistic corrections have been neglected in eq. (5). They can be included perturba- 
tively by inserting higher derivative one-body operators and corrections to the energy-momentum 
relation for external legs. 
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channel where a = —23.714 ± 0.013 fm and tq = 2.73 ± 0.03 fm. The effective range tq and 
the higher coefficients in the momentum expansion are bounded by the range of the nuclear 
potential. 

The coefficients C 2n (/i) are determined by a,r , . . . etc. Their p dependence cancels the 
explicit p dependence in eq. (|) so that the physical scattering amplitude A is p independent. 
The p independence of A implies a set of renormalization group equations for the C*2n's. 



Ai^CW = p2n , (9) 



where, for example, 



Po = ^~C 2 p and P2=¥-C C 2f i, (10) 

47T Z7T 

are the exact beta-functions for C and C 2 . Solving these equations with boundary conditions 
supplied by eqs. (|7|||) yields 

^ = S fc^) - d G ^ = l(qr^) 2 '-»- (») 

In a system with a short range interaction characterized by the momentum scale A, the 
natural size for the parameters a, ro, . . . is set by the scale A. When a is of natural size (i.e. in 
magnitude of order 1/A), a simple power counting like that used in chiral perturbation theory 
for pion self interactions is possible. It is then convenient to take p = which corresponds 
to the MS scheme. The amplitude A can be expanded in powers of p/A with each loop 
contributing a factor of (Mp/4ir) and each vertex a factor of, C2 n p 2n ~ (4np 2n )/(MA 2n+1 ). 
Consequently the leading contribution of order p° comes from Co at tree level. At order 
p 1 a one loop diagram with two Co's contributes. At order p 2 one must include a two loop 
diagram with Co at each vertex and a tree level diagram with a vertex proportional to C 2 , 
etc. 

In contrast, for the special case where \a\ is very large and p = 0, the coefficients are very 
large, C 2n ~ (47ra n+1 )/(MA n ), and the momentum expansion used for \a\ ~ 1/A is valid 
only over a small momentum range. This difficulty is cured by taking p nonzero in which 
case the coefficients Cm need not be large. As \ap\ — > oo one finds, C 2n ~ 47r/(MA n /x n+1 ). 
If one were to set p = A, this would be similar to the result of [11.] where a momentum cut 
off A is used. 

The p dependence induced by the additional finite subtraction in the PDS scheme makes 
it possible to establish a power counting that allows a systematic expansion of the scattering 
amplitude over a much larger range of momentum. Choosing p ~ p, all Feynman diagrams 
with Co's at each vertex give a contribution of the same order. This occurs because Co is 
of order p~ x while each loop integration is of order p. With p ~ p the coefficients Cjn are 
of order p~( n+1 ) and C 2 „p 2n ~ p n ~ l . Consequently only C must be treated to all orders. 
The other coefficients can be treated perturbatively. The scattering amplitude thus has an 
expansion of the form 

A = A-i + A Q + Ai + ... (12) 
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Fig 1. The Feynman graphs contributing to A-\ and A®. Each of the resummed vertices 
is a factor of Co, while the vertex marked with a square is a factor of C2P 2 ■ Summing the 



graphs yield the expressions in eq. ^Tjj) 



where A q is of order p q . The leading amplitude is given by the sum of loop graphs with 
Co vertices, while subleading contributions correspond to perturbative insertions of the mo- 
mentum dependent interactions, dressed to all orders by Co- The first two terms in the 
expansion of A are 



A. 



-C 



and Ar 



-C 2 p 2 



2 > 



(13) 



which arise from the graphs shown in Fig. 1. 

As we are interested in momenta p ~ m, explicit pion fields must be included in the 
effective Lagrangian. The pion couplings are determined from the usual chiral Lagrangian for 
pion-nucleon interactions. Since we are considering nucleon-nucleon interactions in a theory 
where the nucleons are treated nonrelativistically, it is convenient to consider separately 
the pion interaction which gives rise to a static NN potential and the radiation pion field 
responsible for the remaining effects. This is similar to the treatment of the photon field in 
nonrelativistic QED [10,18]. The resulting effective field theory amounts to expanding the 
dependence of nucleon propagators on the 3-momentum of the pion fields when the energy 
loop integrations get a contribution from the pole in the pion propagator. In the 1 S'o channel 
the momentum space potential from one pion exchange is 



K(p,p) 



9a 



ml 



2f 2 W + ml 



(14) 



where q = p' — p, gA = 1-25 and / = 132 MeV is the pion decay constant. It gives an order 



p , spatially extended, nucleon-nucleon interaction since we are treating p 



~ m m 



Hence 



including pions does not change the order A-\ term, and the graphs summed at leading 
order are still those in Fig. 1. For the order p° contribution there are additional terms from 
pion exchange: aside from the contribution to Aq shown in Fig. 1, one must also include the 
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Fig 2. The five sub-graphs contributing to Aq computed in eq. ftUfy. The gray blob is defined 
in Fig. 1, while the dashed line is the exchange of a potential pion. 



exchange of one potential pion dressed by Co to all orders (radiation pions do not contribute 
at this order). 

The expression for Aq is conveniently expressed as the sum of the five graphs shown in 
Fig. 2; Aq is the contribution to Aq in eq. (|13|), while 



(II) 



A 



A UI) 



S\q 



A (v) 




Ap 2 




g\ / m^MA-i 
P \ 

9a f m-nMA 
2f 2 \ 

DoTTll 



tan 



i tan 




(15) 



Note that .Ao has a logarithmic dependence on fi. Since m 2 is proportional to the light 
quark masses we are required to include the vertex in A^ proportional to -D 2 m^ 0. We 
have absorbed into A^ some of the finite, /i-independent, part of Aq that arises in PDS 
(e.g. the part involving Euler's constant and a logarithm of 4tt). In our power counting p, fi 
and m n are considered to be of the same order. Comparing A-\ with Ao indicates that the 
expansion parameter is p/A^N where 



A 



NN 



{8nf/g 2 A M) ~ 300 MeV. 



(16) 



3 The presence of such a divergence Q indicates that the power counting suggested by Weinberg 
is not consistent. The influence of this operator is not subdominant to one pion exchange when 
Cq is treated to all orders. Its presence is needed to get a subtraction point independent amplitude. 
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At leading order in our expansion, for large scattering length, Co(m n ) = —Air/(Mm n ) = 
—3.7 fm 2 . Our power counting leads us to expect that \C 2 (m^)\ and \D 2 (m n )\ are of order 
An/(MA NN ml) ~ 3.5 fm 4 . 

The inclusion of the pions changes the renormalization group equations for Co and C 2 
from what is given in eq. ([y]). To the order we are working the new beta-functions are 

M f o C g 2 A \ . _ M . 

Po = ( C{ ) + -^J H and /? 2 = — C C 2 /i. (17) 

In addition the constant D 2 (fi) satisfies the renormalization group equation 

* D , = !L IW , lt+ us. 

The solutions to these renormalization group equations are the running coupling constants 
Co(/x), C 2 (/j) and D 2 {ji). When expressed in terms of the running couplings, the full ampli- 
tude A is independent of /x; however, as in perturbative QCD, our expression for A-i + Ao 
will be \i independent only up to the order we are working, p°, with residual /i dependence 
canceled by higher order terms in the expansion of A. The free parameters can be taken to 
be the value of the running coupling constants at a particular scale, which we choose to be 
\i = m n . 

All that remains is to fit our free parameters Co(m„-), C2(m n ) and D 2 (m 7r ) to data. We 
compute the phase shift S from the exact expression 

5=iln(l + A0 , (19) 



2i \ 2tt 

expanding both sides to a given order with 5 = \ 5^) + ^(i) + •••)• For example, to subleading 
order, 

We could determine the free parameters analytically by requiring that the theory exactly 
reproduce effective range theory at low momentum. If we do this, our result reproduces the 
data very well up to center of mass momentum p ~ 150 MeV. However, better agreement 
with the data at large momenta results from a fit to the phase shift from the Nijmegen 
partial wave analysis |17j optimized over the momentum range p < 200 MeV. The result of 
this procedure are the values 

C (m n ) = -3.34 fm 2 , D 2 (m n ) = -0.42 fm 4 , C 2 {m n ) = 3.24 fm 4 . (21) 

and the phase shift plotted in Fig. 3. (The individual values of Co and D 2 are not particularly 
meaningful since up to terms higher order in our expansion the amplitude can be written in 
terms of the linear combination Co + -D 2 m 2 .) As is apparent from Fig. 3, the agreement of 
the phase shift with data is excellent at quite large values of p. Furthermore, the coupling 
Co(^Tr) is close to its leading order value (in the limit of large scattering length), —3.7 fm 2 
and C 2 (m n ) is at the expected size, suggesting that our expansion is valid in this channel. 
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Fig 3. The phase shifts 5 in the x Sq channel, and 5q, E\, and 62 in the 3 5*i — Z D\ coupled 
channels, plotted in degrees versus center of mass momentum p. The dot-dashed line rep- 
resents the leading p~ l order calculation for 5 and 6q; at this order one finds S\ = 62 = 0. 
The dashed lines are the results from the order p° calculation. Solid lines are results from 
the Nijmegen partial wave analysis of scattering data J7^/ . The fits involve three parameters 
in the x 5o channel and three parameters in the 3 Si — 3 D 1 coupled channels. 



However, for p > 100 MeV the magnitude of the ratio Aq/A—x is greater than ~ 0.5 and it 
is difficult to justify the approximations we have made, e.g. neglecting terms suppressed by 
(A/A-1) 2 . 

We have performed a similar analysis in the coupled 3 S\ — 3 Di NN scattering channels. 
In this channel the amplitude A is a 2 x 2 matrix with elements Al,v- The S-matrix in this 
channel is usually expressed in terms of two phase shifts, 5q, 62 and a mixing angle ex, 

.pM_ . _ ( e 2iSo cos 2^i - 1 ze i(5o+<52) sin2£i \ , , 

2tt [ ie^+fc) sin 2e 1 e 2ife cos 2e x - 1 ) { ] 

The scattering length is large in the 3 5i channel (a = 5.423 ± 0.005 fm) and the power 
counting is analogous to that of the 1 Sq channel. At leading order, p^ 1 , A is expressed 
in terms of a single parameter that is a linear combination of coefficients of four-nucleon 
operators with no derivatives. Thus A02 = A20 = A22 — 0, which implies that e± = 62 = 0. 
At next order, p°, there are two new parameters analogous to C2 and D 2 that occur in the 



S 



expression for Aqq. In a future publication |I| analytic formulae for the elements All' 
will be given to order p°. Results of fitting the three parameters that occur in Aqq to the 



Nijmegen partial wave analysis []T7 are presented in Fig. 3. Note that we have worked to 
subleading order in calculating 8 , but only leading order in S\ and <5 2 , which explains the 
relative quality of the fits. After fitting 8q there are no free parameters in the predictions 
for E\ and 62 at this order. 

In conclusion, we have shown how to describe NN scattering in a dimensionally regulated 
theory with a consistent expansion, and have demonstrated the failure of Weinberg's power 
counting for these systems. The subtraction scheme we have introduced is useful for treating 
finely tuned theories with power-law divergences such as the standard model with a light 
Higgs, and the Nambu-Jona-Lasinio model, as well as any nonrelativistic theory with a 
large scattering length arising from a short-range interaction. In a future publication we 
will present other applications of the methods developed here, including a discussion of 
other partial waves, inelastic processes, properties of the deuteron, and the iV-body problem 
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